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ABSTRACT 

We present the results of K-band imaging polarimetry of seven 3CR radio galaxies 
with 0.7 < z < 1.2. We find strong evidence for polarization in three sources: 3C 22, 
3C 41 and 3C 114. Of these, 3C 41 shows strong evidence of having a quasar core 
whose infrared light is scattered by dust. We also find some evidence for polarization 
in 3C 54 and in 3C 356. The two pointlike sources (3C 22 and 3C 41) and the barely- 
elongated 3C 54 appear to have of order ten per cent of their K-band flux contributed 
by scattered light from the active nucleus. We conclude that scattered nuclear light 
can form a significant component of the near-infrared light emitted by high-redshift 
radio galaxies, and discuss models in which the scattering particles are electrons and 
dust-grains. 
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1 INTRODUCTION 

A popular model for high -redsh ift radio galaxies was pro- 
posed by Lilly & Longair (1984) a decade ago. Their study 



of the K-z Hubble diagram for 3CR radio galaxies showed a 
tight correlation: the K-magnitudes exhibit a constant dis- 
persion up to z J; 1, but there is evolution with redshift 
such that galaxies at z ~ 1 are about 1 mag more lumi- 
nous than at z = 0. Such a result is explicable if the ob- 
served K-band light originates in stars, and radio galaxies 
undergo a burst of star formation at a very early epoch, fol- 
lowed by passive evolution. A similar model has been sug- 
gested for radio-q uiet elliptical galaxies (Eggen, Lynden-Bell 
fc Sandage 1962| ). 

This model for radio galaxies was challenged by two 
discoveries, which suggested that a significant part of the 
the optical light originated in the active nucleus rather than 
the host galaxy. It was found that optical light from high- z 



radio galaxies is usually aligned with the radio axis (Cham- 
bers, Miley & van Breugel 1987|; |McCarthy et al. 1987f ~ 



the so-called 'alignment effect' - and that the optical light 
is often polarized with its .E-vector oriented pe r pendicular 
to the radio axis (di Serego Alighieri et al. 198S; Jannuzi & 



Elston 1991 



Tadhunter et al. 1992). The polarization of the 



optical light suggests that it is emitted in a restricted range 
of directions by the active nucleus, and has been scattered 
into our line of sight, by either dust or electrons. 

Lilly & Longair's ( 1984 ) model may still be correct if 
the optical light originates in the active core while the (ob- 
server's frame) K-band light comes from stars in the host 
galaxy. One test of this is whether the alignment effect ex- 



tends to the near-infrared; Dunlop & Peacock (1993) find 
a clear infrared/radio alignment effect for a sample of 3CR 
radio galaxie s - co nsistent with the smaller 3C sample of 



Rigler et al. (1992), which suggested that an infrared align- 



ment effect was p resent , but weak. 

Eales et al. (1997) compared the infrared luminosity 
of 3CR galaxies with that of a sample of B2/6C galaxies 
(selected for their lower radio intensities in approximately 
the same frequency band as the 3CR galaxies). They found 
that B2/6C galaxies are fainter than 3C galaxies at simi- 
lar redshifts, by 0.6 mags in K-band at z ~ 1; their sample 
also confirms Dunlop & Peacock's ( 1993] ) finding that the 
infrared alignment effect is str onges t in sources with high 
radio luminosity. Eales et al. (1997) argue, therefore, that 



part of the infrared luminosity of radio galaxies arises from 
the active nucleus, and so a selection effect will produce the 
strong alignment effect and brighter K-band magnitudes ob- 
served in the most radio-luminous galaxies, viz. 3C galaxies 
at high z. One possible explanation of this result is that a 
significant fraction of a radio galaxy's infrared light emerges 
from the active nucleus in a restricted cone, and enters our 
line of sight after scattering by dust or electrons. There is al- 
ready other evidence that light from the active nucleus, if not 
necessarily scattered, can dominate the K -band light from 
high-z radio galaxies: Dunlop & Peacock (1993) noted that 
radio galaxy 3C 22 appeared to be exceptionally roun d and 
bright, and spectroscopic studies by Rawlings et al. (199E) 
confirmed that in the K-band, 3C 22 has properties more 
similar to those of a quasar than of a radio galaxy. 

We have started a programme of near-infrared po- 
larimetry for two reasons. Detection of infrared polariza- 
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Table 1. Rcdshifts and integration times for our sample of radio 
galaxies. 



Source IAU form z A r (/im) t; nt (min) 

3C 22 0048+509 0.936 1.14 72 

3C 41 0123+329 0.794 1.23 108 

3C 54 0152+435 0.827 1.20 135 



3C 651 U22U+3yV TT7E 032 T2 

3C 114 0417+177 e**5 *89 

3C 356 1723+510 1.079 1.06 99 

3C 441 2203+292 0.707 1.29 144 

Key: z: redshift; A r : rest-frame equivalent of observed-frame 
2.2/im; i; llt : total integration time summed over all waveplate 
settings. 

tion would reinforce the evidence that scattered non-stellar 
light makes a significant contribution to the infrared in ra- 
dio galaxies. Also, K-band polarimetry provides a valuable 
complement to optical polarimetry: a long wavelength base- 
line allows a more powerful test of whether the scattering 
centres are electrons, with a wavelength-independent scat- 
tering cross sectio n, or dust, which p referentially scatters 
short wavelengths (Gimatti et al. 1993). 



2 OBSERVATIONS AND REDUCTION 
PROCEDURE 

2.1 Our sample 

We observed a sample of seven high-redshift 3C radio galax- 
ies at redshifts 0.7 < z < 1.2. These form a representative 
sample of the different morphologies present in this redshift 
band: 3C 22 and 3C 41 are bright and appear pointlike; 3C 
114 and 3C 356 display complex knotted morphologies with 
large scale alignments between the K-band structure and the 
radio axis; 3C 54, 3C 65 and 3C 441 are faint sources with 
some indication of K-band structure. Of these three faint 
sources, 3C 54 displays an align ment between the K-ban d 
morphology and the radio axis (Dunlop & Peacock 1993), 



3C 65 shows no prefe rred direction in its H-band structure 
(Rigler & Lilly 1994), and 3C 441 has a broad-band opti- 



cal polariz ation which is roughly perpendicular to its radio 
structure (Padhunter et al. 1992). Pable jj] lists the sources 
observed with their redshifts and the rest-frame wavelength 
of the observed light. 



2.2 Instrumentation 

We had the first run after the commissioning run of a new 
instrument, IRPOL2, installed at UKIRP (the United King- 
dom Infrared Pelescope, Hawaii). Phe IRPOL2 polarime- 
ter consists of a rotatable half-wave plate and a Wollaston 
prism, working in conjunction with the IRCAM3 InSb ar- 
ray detector. We u sed IRCAM 3 at the default pixel scale, 
0.286 arcsec/pixel ([Aspin 1994]), with a K-band filter. Phe 



Wollaston prism causes each source in its field of view to 
appear as a pair of superimposed images with orthogonal 
polarizations, separated by —0.93 pixels in r ight ascension 
and +69.08 pixels in declination (Aspin 1995). 

The design of the instrument is such that when the 



waveplate is set to its 0° reference position, an object to- 
tally linearly polarised with its electric vector at 83° (i.e. 
celestially East of North) would appear only in the upper 
(Northern) image, and an object totally linearly polarised 
at -7° would only appear in the lower image. There are four 
standard offset positions for the waveplate: 0°, 22.5°, 45° 
and 67.5°. Phe IRPOL2 system has negligible instrumental 



polarization (Hough, private communication; Chrysostomou 



199e). We took data on the nights of 1995 August 25, 26, 



and 27. 



2.3 Observing procedure 

For each target object, we took a 'mosaic' of nine images 
with the waveplate at the 0° offset. One image consisted of 
a 60 second exposure (the sum of six ten-second co-adds), 
and the mosaic was built up by taking one image with the 
target close to the centre, and eight images with the frame 
systematically offset from the first by ±28 pixels (8 arcsec) 
horizontally and/or vertically. Phe same source was then 
similarly observed with the waveplate at the 22.5°, 45° and 
67.5° offsets, completing one cycle of observations; hence one 
such cycle took 36 minutes of integration time. Between two 
and five observation cycles were performed over the three 
nights for each target; the total integration time for each 
target is given in Pable [lj Phe times quoted are not exact 
multiples of 36, as in some cases, mosaics were corrupted by 
problems with the windblind, and excluded from our analy- 
sis. An example of a mosaic, 3C 54 and its surrounding field, 
observed with the waveplate at 22.5°, is shown in Figure |l[ 



2.4 Data reduction 

We reduced the data by marking bad pixels, subtracting 
dark frames, and flat-fielding. Flat-fielding frames were ob- 
tained by median-combining the nine images of each mosaic 
without aligning them, and normalizing the resulting image 
by its mean pixel value. In order to align each set of nine 
images, we chose a star present on each frame, and mea- 
sured its position with the APPHOT. center routine in the 
iraf package. Using these positions, the nine images were 
melded into a single mosaic image. 

We performed photometry on the pairs of images of 
field stars and of radio galaxies in each mosaic using AP- 
PHOT. phot from the iraf package. For each target, we chose 
one mosaic arbitrarily to determine the best aperture. Using 
this mosaic, we performed photometry on the two images of 
the source at a series of radii increasing in unit pixel steps. 
Phe measured magnitude in each aperture decreases as the 
light included in the aperture increases; we found the first 
pair (r, r + 1) of radii where the change in magnitude was 
less than half the measured error on the magnitude. We then 
used the photometry at the next radius, (r + 2). In this way, 
we hope to include all the source light but as little as possible 
of the surrounding sky. Where the chosen apertures differed 
for the two images of the source, we adopted the larger. The 
chosen aperture was then used on all the mosaics containing 
the target to obtain a set of intensities at the corresponding 
waveplate positions. 

The APPHOT.PHOT routine corrects for the sky bright- 
ness by measuring the modal light intensity in an annulus 
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Figure 1. K-band image of 3C 54 (at either end of the bar) and 
surrounding field. (North is at the top, East at the left. The image 
was composed by mosaicing a series of nine 60-second exposures 
taken with the waveplate offset at 22.5°.) 



around the target; the position of the annulus was chosen in 
each case such that the outer radius did not extend to the 
nearest neighbouring object, and the inner radius was nor- 
mally set one pixel greater than the photometry aperture. 
(Where we attempted to perform photometry on a knot 
within a larger structure, the inner radius of the annulus 
was set sufficiently large to exclude all the knots comprising 
the object.) This procedure results in pairs of photometric 
intensities for each target. 

2.5 Confidence intervals for polarization 

Using the photometric intensities from all the observing 
cycles, we pooled the data and calculated the intensity- 
normalised Stokes Para meter s, q and u, following the 
method of Clarke et al. (1983). The q and u (as quoted in 
Table ^) are in the instrumental reference frame, such that 
a positive q and zero u would indicate polarization at 83° 
East of North: q is defined by subtracting the lower beam 
intensity from the upper with the waveplate at the 0° refer- 



ence position, and u similarly with the waveplate at 22.5°. 
The Stokes Parameters were t hen used to derive a measured 
degree of polarization, p — y^q 2 + u 2 . 

Since p is positive in the presence of noise even for an 
unpolarised source, there has been some debate about the 
best metho d of obtaining an 'unbiase d ' estimate of t he po- 
larization ( Simmons fc Stewart 1985 ; Leyshon 1997 ). Fol- 
lowing best practice, we assume that the Stokes Parameters 
q and u are normally distributed, and that the errors on the 
two parameters are comparable. Under these assumptions, 
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Underlying polarization, a 



Figure 2. Illustration of contours in the a — m plane. See text. 



it can be shown (Wardle & Kronberg 1974; Vinokur 196 
that, if the underlying polarization of a source is po and the 
error on a measurement is a, then the probability distribu- 
tion function for the measured polarization values p will be 
a Rice distribution. Defining normalised variables a =po/a 
and m — p/er, we can write this distribution as 



F(m, a) = m. exp 



Jo(ma) 



(1) 



where m ^ and Jo is the modified Bessel function. It fol- 
lows, trivially, that the integral 



m— mc (a) 



F(m, a). dm = C 



(2) 



for ^ C ^ 1 represents a confidence interval (0, mc) such 
that there is a chance C of a measurement of the polarization 
lying within these bounds. 

We need to invert the probability distribution to es- 
timate a given our meas ured m. Following the method of 
Mood, Graybill & Boes ( 1974 ) we can solve Equation |^ for 
mc for given values of C and a, and so plot contours of con- 
stant C on the a-m plane, as illustrated in Figure |^. Clearly 
two contours C\ , C2 allow us to read off a C2 — Ci confidence 
interval (mt,,mxj) for given ao- Mood et al. ( ]l974[ , Ch. VIII, 
§4.2) show that the same contours also provide a C2 — C\ 
confidence interval for the underlying polarization (aL,au) 
for a given measurement mo- 

As well as providing a general interval (ol, ou), the con- 
tours can be used to determine the probability that the un- 
derlying polarization is non-zero. The probability that a — 
exactly is finite rather than infinitesimal, since some sources 
may be unpolarised. There will be some contour C m , which 
cuts the m-axis at our measured mo, and so divides the do- 
main of a into a > with confidence C m , and a = with 
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confidence 1— C m . This confidence can be found analyticaliy, 
since Equation ^ can be integrated at a — 0; and it follows 
that the probability that a > 0, given some measured value 



C m = 1 - exp( 



-m 2 /2). 



(3) 



Equation g can be used to determine the probability 
that a source is truly polarised, and the probabilities so de- 
rived for each source are listed in Table |^. 

For a given measurement mo, it is also possible to use 
appropriate pairs of contours to find a set of 67 per cent (la) 
or 95 per cent (2a) confidence intervals for our best estimate, 
a. For any given required confidence, there will be a continu- 
ous set of pairs of cont ours which could be chosen; following 
Simmons & Stewart (1985) we choose the unique interval 



whose width is minimized by the additional constraint 

F(m L , a) = F(mv, a). 

The errors on the polarizations quoted in Table |^ are the la 
errors derived by this method. 



(4) 



2.6 Point estimates for polarization 

To obtain a point estimate of t he polarization, we again 
follow Simmons & Stewart (1985), who have tested various 
estimators &x for bias. They find that when a & 0.7, the 
best estimator is the 'Maximum Likelihood Estimator', 5ml; 
when a <; 0.7, the best estimator, fiwK, is that traditi onall y 
used by radio astronomers, e.g. Wardle fc Kronberg (1974) 



It ca n be shown that «ml < bwk Vm (Simmons & Stewart 



1985) 



We do not know, a priori, whether a will be greater or 
less than 0.7, and hence which estimator to use. We distin- 
guish three cases: (i) both estimators indicate fix S5> 0.7, so 
we must use 5ml; (ii) both estimators indicate &x 0.7, so 
we must use aw«; or (iii) in the region of ax ~ 0.7, each 
estimator will indicate that it should be preferred over the 
other. We therefore wish to construct an estimator a which 
is continuous with m and takes these cases into account. 
This can be done using the inequality oml < Owk Vm. 

Case (i): Let mwKmin be the measured value of polar- 
ization which would cause an estimate 5ml < Q-wk = 0.6. 
Clearly if m < mwKmin then a < 0.7 and the Maximum 
Likelihood estimator is cert ainly the most appropriate. We 
find that mwKmin = 1.0982 (Leyshon 1997) and in this case, 
the Maximum Likelihood estimator will in fact be zero. 

Case (ii): Let mMLmax be the measured value of polar- 
ization which would cause an es timat e 0.8 = <2ml < &wk- In 
this case, Wardle & Kronberg's (1974) estimator will clearly 
be the most a ppropriate. We find that m-MLmax = 1.5347 
( |Leyshon 1997] ). 

Case (iii): Between these two extremes, if our measured 
polarization satisfies mwKmin < m < niMLmax, we can inter- 
polate such that our best estimate is 



rn — niwKn 



WlMLrr 



•CtML ' 



TlMLmax — myVKrr 



■%K- 



(5) 



Thus we can obtain both a 'debiased' point estimate 
and a confidence interval for a, and hence for the degree of 
polarization aa. 



2.7 Estimating the angle of polarization 

The position angle, <f> = 83° +(f> m , of the E- vector is obtained 
by calculating 



2<f>r, 



tan 1 (u/q), 



(6) 



and noting the appropriate quadrant from the signs of q 
and u. There are two methods of obtaining the error on 
this angle. We can propagate the errors on q and u through 
Equation [| yielding an asymmetric error; or the symmetric 
67 per cent confidence interval may be constructed, using 
the best estimate of a to calculate t he relevant dis tribution 
function ( Wardle fc Kronberg 1974 |Vinokur 1965 ). We cal- 
culated error bars by both methods, and selected the greater 
error in each case to quote in Table |^. 

2.8 Extinction corrections 

The interstellar medium may impose some linear polariza- 
tion on light transmitted through it. Empirically, the polar- 
ization p depends on wavelength, such that 



P/Pmax = exp[— 1.7 A n 



: In 2 (A n 



c/A)] 



(7) 



where A max is the wavelength (explicitly in micrometres) 
at which the polarization peaks, us ually around 0.5pm, and 



empirically i n the range 0.3-0.8 p m (Serkowski, Matthewson 
Wilking et al. 198C| ) . The relationship is found 



fc Ford 1971 



t o hold up to around 2um , and is adequate for our purposes 
flMartin fc Whittet 199c| ). 

In general, p max for a given set of galactic co-ordinates 
is not known. But suppose we take the ratio of polarizations 
in two wavebands, V and K, and rearrange: 



Pk = Pv exp 



-3.4A n 



vAkAv 



(8) 



but for 0.3 < 



Hence pk = c.pv where c depends on A 
A max < 0.8 we find 0.15 < c < 0.30. 

An empirical upper limit for py (e xpressed as a percent- 
age) is given by Schmidt-Kaler ( |l95£ ) as pv = 9Eb-v and 
typically, py 



4.5£b- 



It is well established (Sava ge & 
Mathis 1979| ; |Koorneef 198^ ; |Rieke fc Lebofsky 1985] ) that 
the ratio of total to selective extinction is Av/Eb-v ~ 3; 



and so we can use the extinctions Ab (Burstein fc Heiles 



1982, figures obtained from the NED database) to obtain 
Eb-v =Ab-A v = Ab/4. 

Taking the maximum values, c = 0.3 and pv = 9Eb-v, 
we find an upper limit for infrared polarization pk = 0.7j4b- 
The NED values for Ab and the corresponding upper limits 
on pk are given in Table ^[ 



3 RESULTS 

3.1 Have we detected polarization? 

Equation |^ allows us to quantify the probability that a given 
object is polarised. The probabilities of each object being 
polarised are listed in Table ^. Three of our seven sources 
have a 95 per cent or better probability of being polarised; 
and of these, 3C 22 and 3C 41 are polarised at the 3 per 
cent level, and 3C 114 at the 12 per cent level. 

The number of prominent starlike objects (in addition 
to the target) featuring on the mosaics varies between one 
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Table 2. Total B-band extinctions and upper limits on K-band 
interstellar polarization for our sample. 



Source 


Ab 


PK,max 


3C 22 


1.09 


0.76 


3C 41 


0.17 


0.12 


3C 54 


0.37 


0.26 


3C 65 


0.16 


0.11 


3C 114 


1.26 


0.88 


3C 356 


0.10 


0.07 


3C 441 


0.34 


0.24 



Key: Ab: Blue-band ex 
from Burstein & Hcilcs (1198 



on (magnitudes), from NED, derived 



; Pk, 



: maximum interstellar con- 



tribution to K-band polarization (per cent) as described in the 
text. 

and seven, depending on the target. Where possible, we have 
performed polarimetry on the stars; out of the 21 stars so 



observeil, uiily uue lias a giealm Lhaii 95 pel cent piubabiliL i y 
of bein g polarised. Thia object waa a bright atarlikc object 

on the mosaic containing 3C 114, but is only polarised at 
the 0.7 per cent level, which is explicable by the interstellar 
medium (see Table |^) . 

Within the bin of sources having a probability 80-95 per 
cent of being polarized, fall three further stars; of these, one 
is extremely faint, and another appears to be polarised at 
only the 0.3 per cent level. The third falls on the same mosaic 
as 3C 54, and appears to be polarised at the 5.6 ± 2.6 per 
cent level, with a 94 per cent chance of the polarization being 
genuine. This star, however, straddles the edge of three of 
the nine component frames of the mosaics, so the validity of 
the result is called into question. Two of our sources also fall 
in the 80-95 per cent probability bin: 3C 54 itself, polarised 
at the 6 per cent level, and 3C 356, at the 9 per cent level. 

Given that 17 out of 21 stars, but only 2 out of 7 sources, 
have a probability of less than 80 per cent of being polarised 
at all, we feel confident of having detected polarization in 
three 3C sources, and possibly in a further two. 



3.2 Individual objects 

3C 22 

3C 22 has a 95 per cent chance of truly being polarised, and 
the debiased polarization is 3.3 ± 1.4 per cent. We expect 
that no more than 0.8 per cent is due to the interstellar 
medium; most of the polarization is therefore intrinsic to 
the source. The orientation of the B-vector is +27^7!° East 
of North. Our K-band image shows no evidence for extended 
structure. 

The radio posit ion angle of 3C 22 is 102 
Kapahi 



( Schilizzi 



Neff 1982). The errors on the measured polariza- 
tion orientation are large, but the figures suggest that the 
true direction is more likely to be perpendicular to the radio 
axis, than parallel to it. As we have alr eady seen, 3C 22 is 
suspected of being an obscured quasar (Dunlop & Peacock 



199S 



; Lon- 



awlings et al. 1995 ); according to data in Lilly , 
gair ( P-984T ), its K-band magnitude (15.67±0.10) its brighter 
than th e mean K-z relationship by about 0.9 mag. A perpen- 



dicular alignment of the radio axis and iJ-vector orientation 
would be expected for light from the active nucleus travel- 



ling along the radio axis and then scattered towards us by 
either dust or electrons . Jannuzi (private communication) 
(Elston & Jannuzi 1997) has performed imaging polarime- 
try on 3C22 at shorter wavelengths, and reports 3<j upper 
limits in V and H of 5 per cent and 3 per cent respectively. 



3C 41 

3C 41 has a 98 per cent chance of having an underlying 
polarization, which we measure to be 3.1 ± 1.1 per cent. 
Our upper limit for extinction-induced polarization is only 
0.1 per cent, so we are confident of having detected intrinsic 
polarization in this object. The orientation of the B-vector 
is +45^14° East of North. Our K-band image shows no evi- 
dence for extended structure. Like 3C 22, Lilly & Longair's 
( |l984| ) data shows that, at K = 15.95 ±0.10, 3C 41 is signif- 
icantly brighter than the mean K-z relationship, by about 
0.6 mag. 



Jannuzi (private communication) (Elston & Jannuzi 



1997) have firm V and H band polarizations for this source: 
at V, the polarization is 9.3 ± 2.3 per cent at +58°± 7°; 
at H, the polarization is 6.6 ± 1.6 per cent at +57°± 7°. 
The .E-vector orientations in the three wavebands are con- 
sistent with_^ne_another. The radio position angle of 3C 41 is 
147° (Longair 1975), so we have a very good perpendicular 
alignment between the radio and polarization axes. 



3C 54 

The probability that 3C 54 is polarised is 94 per cent; our 
measured value of polarization is 5.9 ± 2.6 per cent at 
+32l 7 g° East of North. Dust is not expected to contribute 
more than 0.3 per cent. We have no independent polarimetry 
for t his object. Th e position angle of its radio structure is 
24° ( Longair 1975 ); our measurement shows that the polar- 
ization orientation is more likely to be parallel to the radio 
axis, than perpendicular to it, b ut th e latter case cannot be 
ruled out. Dunlop & Peacock's ( 1993| ) K-band contour map 
shows a slight extension to the south-west, which was also 
just discernible in our image. This would be roughly parallel 
with the radio axis. 



3C 114 

An image of 3C 114 is shown in Figure ^. It consists of at 
least four knots, where the brightest knot forms the knee 
of a P-shaped structure. Only the knee knot proved bright 
enough to analyze on its own; Table ^ includes results for 
both that knot and the structure as whole. There is a prob- 
ability in excess of 99 per cent that there is genuine polar- 
ization in both the knee knot and the overall structure. 

We measured polarization in the knee of 5.1 ± 1.7 per 
cent, at +llllg3° East of North. Overall, the whole object 
has a polarization of 11.7 ± 3.0 per cent at +73^20° • The 
extinction contribution could be as high as 0.9 per cent, but 
our detections of polarization are much higher than this, so 
the polarization appears to be intrinsic. 



The radio position angle of 3C 114 is 54° (Strom et al 



199C). Our measurements of the i5-vector orientation sug- 
gest that an alignment parallel to the radio and optical axes 
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Table 3. Observational results from polarimetry of 7 radio galaxies in K-band. 



Source 




)(%) 


«K)(%) 


r(") 


P±a P (%) 


2cr.UL 


Prob 


en 




3C 22 


-1.3 


(1.4) 


-3.2 (1.4) 


2.6 


3.27 ± 1.38 




0.95 


27.0 


-71.4 +16.9 


3C 41 


+0.8 


(1.1) 


-3.1 (1.1) 


2.3 


3.09 ± 1.11 




0.98 


45.4 


-14.3 +78.6 


30 54 


-1.4 


(2.5) 


-6.0 (2.5) 


1.0 


5.90 ± 2.55 




0.94 


31.6 


-79.1 +17.3 


3C 65 


-4.3 


(4.2) 


-1.2 (4.0) 


2.9 


2.15 ± L 4.49 


9.7 


0.42 


0.6 


±71.4 


3C 114 


+11.2 


(3.1) 


-4.1 (2.7) 


3.6 


11.74 ± 3.02 




0.99 


72.9 


-20.6 +33.9 


•Knee 


+3.0 


(1.6) 


+4.3 (1.8) 


1.7 


5.11 ± 1.73 




0.99 


110.6 


-63.2 +13.5 


3C 356 


-9.6 


(4.9) 


+3.5 (5.6) 


4.0 


9.04 ± 4.94 


16 


0.85 


172.0 


±22.1 


•North 


-13.3 


(7.9) 


+4.6 (8.4) 


2.3 


13.02 ± 7.99 


41 


0.62 


163.5 


±25.0 


•SE 


-10.1 


(9.4) 


-18.9 (16.7) 


2.6 


18.99 ± 15.31 


24 


0.78 


23.9 


-54.2 +33.0 


3C 441 


-1.7 


(2.3) 


-0.8 (2.3) 


3.4 


1.04 ± L 2.35 


4.7 


0.26 


6.2 


±77.3 



Key: q ± a q , u ± a u : normalized Stokes Parameters (per cent) with respect to 83° E of N; r: radius of photometry aperture 
(arcseconds) ; P±<rp: percentage polarization (debiased) with la error — the la lower limit for polarization is zero); Prob: 
the probability that there is underlying polarization, given by Equation pt 2ct.UL: 2a upper limit (in per cent) for polarization 
in objects unlikely to be polarised; 8 ± erg: Electric vector orientation E of N (°). 



4 DISCUSSION 



Figure 3. K-band structure of 3C 114. (North is at the top, East 
is on the left.) 



is more likely than a perpendicular alignment, but the errors 
do not permit a more definite conclusion. 



3C 356 

We have no strong evidence for polarization in 3C 356, but 
we note that the structure has three lobes in its K-band 
image, at the North, the South-East and the South- West. 
Only the North and South-East knots were bright enough 
to permit polarimetry, as recorded in Table ra. 



4.1 The fractional contribution of quasar light 

Radio galaxies are thought to have quasar nuclei at their 
cores, but to be oriented such that no direct radiation from 
the core can reach us. Nearby radio galaxies are known to 
have the morphology of giant ellipticals; if this is also true at 
higher redshifts, the stellar populations of our seven sources 
will consist of mostly old, red, stars. The total light received 
from our sources will be a combination of starlight, light 
from the active nucleus scattered into our line of sight, and 
nebular continuum emission. 

We denote by $w, the fraction of the total flux density, 
Fw, in a given waveband, W, which originates in the active 
nucleus. We expect that in the visible wavebands, $u,b,v 
will be a significant fraction of unity. From our observations, 
we wish to determine whether <E>k is small, or whether a 
significant component of the stellar-dominated infrared also 
arises in the active nucleus. 

Where the magnitude, W, has been measured in a given 
waveband, the total flux density can be calculated: 



Fw =Fo(W = 0).10" 



(9) 



if we denote the flux density scattered into our line of sight 
from the quasar core by Fq,w, then $w = i*Q,w/-Pw- 

The optical flux density of quasars can be modelled well 
by a power law of the form CM~ a where a, the 'spectral 
index', is of order unity. Since the efficiency with which a 
given species of scattering centre scatters light may depend 
on wavelength, we denote that efficiency by fw, and the 
scattered quasar component can be expressed 



-Fq.w = C.jw-v 



(10) 



It will be computationally convenient to define an 'un- 
sealed model light ratio', (f>w, as 



4>w 



fw-v 



F (W = 0).10" 



(11) 



then the actual model light ratio will be $w = C<j)w- There 
is clearly an upper limit set on C by the fact that $w niay 
not exceed unity in any waveband; hence C < l/(j>w- To 
allow for errors in the observed magnitudes, AW, the max- 



© 0000 RAS, MNRAS 000, 000-000 



K-band polarimetry of high-redshift radio galaxies 7 



imum permissible value of C in a model can be taken to 
be 



C max = min 



F (W = 0).10 



-0.4(W-AW) 



,vw, 



(12) 



where the value for a and the choice of scattering centre 
(hence of /w) depend on the model. 



4.2 The dilution law for polarization 



Following Manzini & di Serego Alighieri (1996) (hereafter 
MdSA) we assume that only the scattered component of the 
light from radio galaxies is polarised. It can readily be shown 
that when linearly polarised light is mixed with unpolarised 
light, the overall degree of polarization is directly propor- 
tional to the fraction of the total intensity contributed by 
the polarised component. Hence defining Pq,w as the intrin- 
sic polarization produced by the scattering process, and Pw 
as the observed polarization after dilution, it follows that 



$w = Pq,w/Pw = Pw/Pq; 



(13) 



If we know the restrictions on possible values of the 
intrinsic polarization Pq,w, we can use our corresponding 
measurements of Pw to put limits on $w for the measured 
sources. The appropriate restrictions depend on the nature 
of the scattering centres. 



If the scattering centres are electrons (Fabian 1989), 
then Thomson scattering will take place: the effects of the 
geometry and of the wavelength can be treated indepen- 
dently. The spectral energy distribution of the light scattered 
in a given direction is independent of the scattering angle: 
Pq,w = Pq will be a constant. The degree of polarization 
of the scattered light is given simply by 



cos x 



cos z x 



(14) 



where \ i s the scattering angle. For an AGN observed as 
a radio galaxy, we assume an orientation 45° < \ < 90° > 
whence 1/3 < Pq < 1 V W. 

The case where the scattering centres are dust grains 
has been modelled recently (MdSA); the fraction of the 
light scattered by the dust, /w, and the polarization of 
the scattered light, Pq,w, both depend strongly on wave- 
length. The exact relationship depends critically on the size 
distribution of the dust grains, and the amount of extinc- 
tion they introduce; MdSA provide a series of graphs for 
the variation of /w and Pq,w with rest-frame wavelength 
0.1 /im < A r < 1.0 /im, corresponding to many different 
dust grain compositions and size distributions. At the red- 
shifts of the objects in our sample, the light observed in 
the H and J bands originated at rest-frame wavelengths be- 
low 1.0 fim, but the K-band light originated in the region 
1.0 /xm < A r < 1.15 /zm. To accommodate the K-band light 
within our models, we linearly extrapolated MdSA's curves 
to A r = 1.15 /im. 

In models where the smallest dust grains have a radius 
not less than 40 nm, Pq,w approaches zero twice: once at a 
(rest frame) wavelength around 0.2 /im, and again at some 
wavelength between 0.1 and 0.7 /tm which depends critically 
on the dust grain size distribution. But in all cases, Pq,w 
extrapolated into the 1.0 /im < A r < 1.15 /im region gives 



Table 4. Percentage of K-band light estimated to be arising from 
the active nucleus in the 5 polarized galaxies. 



Source 


Pk 






3C 22 


3.3 ± 1.4 


7 ± 4 


8 ± 4 


3C 41 


3.1 ± 1.1 


6 ± 4 


8 ± 3 


3C 54 


5.9 ± 2.6 


12 ± 8 


15 ± 7 


3C 114 


11.7 ± 3.0 


23 ± 13 


29 ± 10 


3C 356 


9 ± 5 


18 ± 13 


23 ± 13 



Key: Pk: measured K-band polarization; &K e : fraction of light 
from quasar according to electron model; <3?Kd : fraction of light 
from quasar according to dust model. Depolarization from mul- 
tiple scattering (both models), and shallower scattering angles 
(dust model only) will tend to increase $k- The polarizations 
used for the compound objects 3C 114 and 3C 356 are those for 
the whole compounds. 

0.3 < Pq.k < 0.5; the intrinsic polarization at V and H is 
lower. 

For the five sources in which we have evidence of K- 
band polarization, we can hence estimate $k under both 
electron and dust models. The values are given in Table ^j. 
For the dust models, we take 1/Pq,w = 2.5 ± 0.5; the error 
takes into account all dust models, and the different redshift 
corrections for the different galaxies, but assumes that the 
scattering angle is 90° . If the scattering angle is less, we as- 
sume that less polarization occurs (see MdSA, Figure 20), 
and hence will be greater than our estimate. For the 
electron models, we multiply the observed polarization by 
1/Pq,w = 2 ± 1; this takes into account all possible orienta- 
tion effects. 

The effects of multiple scattering have been ignored for 
both models; multiple scattering tends to depolarise light, 
and so the true value of under multiple scattering will 
again be greater than our estimate. The only physical mech- 
anism which could cause <&k to be lower than our estimate, 
is polarization of light in transit by selective extinction; and 
as we have already seen (Table ^|), any such contribution to 
the polarization will be small. 



4.3 Modelling the polarization of radio galaxies 

Given a measurement of the magnitude of a radio galaxy, we 
can predict its polarization as a function of wavelength, up 
to a multiplicative constant. Rearranging Equation ll3L and 
employing our 'unsealed model light ratio', we first define 
an 'unsealed model polarization' IIw = Pq,w-0w, and so 
express our modelled polarization as: 



W, modelled 



^w-f&w = C.Pq,w-0w = C.IIw- 



(15) 



To fit a dust scattering model, we can calculate IIw by 
obtaining /w and Pq,w from suitable curves in MdSA. For 
electron models, the wavelength-independent term /w can 
be considered to have been absorbed into the multiplicative 
constant, C, while Pq, also wavelength-independent, can be 
assumed to be its minimum value, 1/3. We cannot sepa- 
rately identify C and Pq; the physical constraints on these 
constants are 1/3 < Pq < 1 and < C < C m ax- If Pq is 
greater than the assumed 1/3, then C will be correspond- 
ingly smaller. 



© 0000 RAS, MNRAS 000, 000-000 



8 G. Leyshon and Stephen A. Eales 



O 



: 

O 
— 

CO 

cd 
O 
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Figure 4. Measured and best-fit model polarizations of 3C 4i 
as a function of rest-frame wavelength. Solid line: dust model, 
a = 1.167; dashed line: electron model, a = 1.733. 



Given a set of N polarization measurements Pw ± ow, 
and a corresponding set of unsealed model polarizations, 
ITw ± ew, based on measured magnitudes, we can calculate 
the deviation of the fit: 



5 = 



N ' (T 

Wt...W N 



c.ik 



2 +(C7.e w ) 2 



(16) 



The best fit is that with the value of C which minimizes 8, 
subject to the physical constraint < C < C ma x. 



A model for 3C 41 

The source for which we had the most data was 3C 41, with 
polarimetry in 3 bands: Pv = 9.3 ± 2.3 per cent; Ph = 6.6 
± 1.6 per cent; Pk = 3.1 ± 1.1 per cent. We attempted to fit 
two models; an electron model and a typical dust model with 
a minimum grain radius of 80 nm. To fit the models to the 
observed polarizations, we tested a discrete series of possible 
spectral indices, —0.5 < a < 2, with a step size of 1/3. For 
each value of a we calculated the 'unsealed polarizations' 
IIv ± ev , IIh ± £h , IIk ± £k ■ We then iteratively determined 
the best fit value of C for each a, and took as our overall 
best fit that combination of a and C which gave the lowest 
6. 

Ma gnitudes for 3C 41 w ere available in 5 bands: J, H 
and K (|Lilly fc Longair 1984j) and the narrow filters g and 
rs (Dickinson, private communication). The H and K val- 



ues yielded direct estimates of the corresponding 'unsealed 
polarizations' H~h and IIk; II v was estimated by linear in- 
terpolation between n g and Il rs . 

Figure ^ shows the measured polarizations for 3C 41 
as triangles (a) and the magnitude-based polarization esti- 
mates, after best-fit scaling, as stars (tV). The lines give the 
error envelope on the modelled polarization (based only on 



the errors on the magnitudes) . The solid lines correspond to 
the dust model, and the dashed lines to the electron model. 

As can be seen from Figure ^j, the model curves lie be- 
low the data point at H, but above those at K and V. The 
shape of the curve depends more strongly on the measured 
magnitudes (and on /w for dust) than on the spectral in- 
dex, and the models for all reasonable values of a will have a 
broadly similar shape; the best fit will necessarily pass below 
the measured point at H, and above that at V. 

Consider the electron model. Figure ^ shows us that 
the theoretical polarization curve for electron scattering is 
concave with respect to the origin, whereas a curve through 
the three data points would be convex; clearly it will not be 
possible to obtain a close fit for the central (H-band) point. 
Fitting the electron model curve, we found that the best fit 
occurred for a — 1.733, with a deviation 8 — 0.430. From 
Table [I| we have $ K = 6±4 per cent for 3C 41. Multiplying 
the polarizations observed at H and V by 1/Pq,w = 2 ± 1, 
we predict 3>h = 13 ± 7 per cent, and <I>v = 19 ± 10 per cent 
in these bands. 

The dust model chosen as typical from MdSA was that 
for a cloud of spherical dust grains, with radii 250 nm > 
a > 80 nm, with the number density per unit dust mass 
following an a~ 3,5 law. This model produced a curve which 
fitted the data points very well. The best fit indicated that 
the optimum spectral index was a = 1.167, for which 8 — 
0.177. 

This dust model was also used to estimate the propor- 
tion of scattered light at shorter wavelengths: $k = 8 ± 4 
per cent, $h = 24 ± 6 per cent, and $v = 155 ± 38 per cent. 
MdSA's curve for polarization as a function of rest-frame 
wavelength predicts a 6 per cent polarization in the observed 
V-band, l ower than the 9 per cent a fter dilution measured b y 
Jannuzi ( private communication] ) (Elston & Jannuzi 1997). 



This is still consistent, within error bars, as long as the true 
value of <E>v for 3C 41 is less than, but very close to, unity; 
the observed V-band corresponds to the near-ultraviolet in 
the rest frame of 3C 41, and it is reasonable (MdSA) to sup- 
pose that the scattered quasar light in that band could form 
in excess of eighty per cent of the total light. 

We noted earlier that the shape of the dust model po- 
larization curve between 0.2 /im and 0.7 /im (rest frame) is 
very sensitive to the choice of dust grain distribution. The 
particular dust model chosen approaches zero polarization at 
a wavelength corresponding to the rs band when redshifted 
into our frame. This causes the 'well' visible in the model 
polarization curve (Figure ^), whose presence is essential for 
the dust model curve to fit the data points closely. 

We also considered other dust models from the selection 
given by MdSA. Of those which differed significantly from 
the 'typical' one considered so far, many of them will not 
predict polarizations in the observed V-band which are suf- 
ficiently high to be reconciled with the observed 9 per cent; 
and those which do, do not possess the deep well needed to 
fit the polarimetry across the spectrum. We conclude, there- 
fore, that the best model for 3C 41 is that of an obscured 
quasar core with a ~ 1.2 beaming its optical radiation into 
a dust cloud, although an electron model cannot be ruled 
out within the error bars. 
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Figure 5. Modelled (o = 1) and measured polarizations of 3C 
22 as a function of rest-frame wavelength. Solid line: dust model; 
dashed line: electron model. 



A model for 3C , 
For 3C 22, we have one firm pol arimetry point (this paper) 



and two upper limits in V and H ( Jannuzi, private communi- 
catio ii[T|( |Elston fc Jannuzi 19971); magnitudes were ava ilable 



in 4 bands including J, H and K ( Lilly fc Longair 1984| 
a crude eye estimate in r ( Riley, Longair fc Gunn 1980| 



and 
The 

measurements and models are shown in Figure ^, with the 
same symbols as Figure |i|; open triangles represent upper 
limits. We have taken a = 1, and normalised the theoretical 
curves to the K-band data point. 

Here the models are inconclusive. In the observed near- 
infrared, both models can easily fit, with some slight scaling, 
within the K-band measurement error bars; both models 
suggest that about 8 per cent of the K-band light arises in 
the active nucleus (Table [|). The error on the observed r- 
band magnitude is so large that both models are consistent 
with the observed V-band upper limit of polarization. 



5 CONCLUSIONS 

We have taken K-band polarimetry for seven 3CR radio 
galaxies, and found a diverse range of results. Out of our 
seven sources, two (3C 65 and 3C 441) display no evidence 
for polarization. For the five sources which do display some 
evidence of polarization, we have estimated the fraction of 
observed K-band light which originated in the active nu- 
cleus. 



The recent finding (Eales et al. 1997) that 3C galaxies 



at z ~ 1 are 1.7 times as bright as the radio- weaker 6C/B2 
galaxies in a similar sample requires that <J?k ii 40 per cent 
(= 7/17) for 3C galaxies, if the scattering of nuclear nonstel- 
lar light is responsible for the brighter K-band magnitudes 
of 3C galaxies. Our results suggests that "!>k is somewhat 
lower. Nevertheless, our results support the hypothesis that 



radio galaxies consist of quasar nuclei embedded in giant 
elliptical galaxies. 

All of the galaxies which appear to be polarised have 
large errors on the orientation of the _E-vector; hence any 
apparent alignment effects are suggestive rather than defini- 
tive. But with this caveat, we note that two sources (3C 54 
and 114) display extended K-band structure, and have high 
polarizations oriented in roughly parallel alignment with the 
radio axis and extension of the optical structure - i.e., in the 
opposite sense to the perpendicular alignment expected un- 
der a simple scattering model. 

Finally, 3C 22 and 3C 41 display significant polariza- 
tion of around 3 per cent, with a polarization alignment 
perpendicular to their radio axes; both appear in K-band 
as pointlike objects. We suggest, therefore, that in these ob- 
jects, infrared light from a quasar core is being scattered 
into our line of sight, and forms a significant part of the 
total K-band flux received from these sources. 

Manzini & di Serego Alighieri (MdSA) comment that 
given the range of possible dust models, 'the wavelength 
dependence of polarization is not necessarily a discriminant 
between electron and dust scattering.' The data available to 
us are insufficient to indicate whether the scattering centres 
in these objects are electrons or dust; it is not possible to 
give an unambiguous fit of the polarization curves with so 
few data points, although 3C 41 does seem to fit a model 
(MdSA) with a minimum dust radius of 80 nm, particularly 
well, and we suggest that it does indeed consist of a quasar 
obscured by dust. 
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